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Abstract — We consider the problem of quantifying the Pareto 
optimal boundary in the achievable rate region over multiple- 
input single-output (MISO) interference channels, where the 
problem boils down to solving a sequence of convex feasibility 
problems after certain transformations. The feasibility problem is 
solved by two new distributed optimal beamforming algorithms, 
where the first one is to parallelize the computation based on the 
method of alternating projections, and the second one is to local- 
ize the computation based on the method of cyclic projections. 
Convergence proofs are established for both algorithms. 

Index Terms — MISO-Interference Channel, Distributed Beam- 
forming, Achievable Rate Region, Pareto Optimal 

I. Introduction 

Traditional wireless mobile systems are designed with the 
cellular architecture, in which neighboring base stations (BSs) 
in different cells try to manage communications for their in- 
tended mobile stations (MSs) over non-overlapping channels. 
The associated inter-cell interferences, from non-neighboring 
cells, are treated as additive background noises. To improve the 
performance of traditional systems, most beyond-3G wireless 
technologies relax the frequency reuse constraint such that 
the whole frequency band becomes available for all cells. As 
such, joint signal processing across neighboring BSs is needed 
to cope with the sttong inter-cell interferences in the future 
cellular systems. 

In this paper, we study a particular type of multi-BS 
cooperation for downlink nansmissions, where we assume a 
scenario with each BS equipped with multiple antennas and 
each MS equipped with a single antenna. Besides, only one 
MS is assumed to be active in each cell at any given time 
(over a particular frequency band). Our problem setup can be 
modeled as a multiple-input single-output (MISO) Gaussian 
interference channel (IC), termed as MISO-IC. 

From an information-theoretic viewpoint, the best achiev- 
able rate region to date for an IC was established by Han 
and Kobayashi in [1], termed as the Han-Kobayashi region, 
which utilizes rate splitting at transmitters, joint decoding 
at receivers, and time sharing among codebooks. The Han- 
Kobayashi region was simplified in (2] and proved to be 
within 1-bit of the capacity region of the Gaussian IC in 0. 
However, in cellular systems, practical constraints often limit 
MSs to only implement single-user detection (SUD) schemes, 
i.e., treating the interference from all other unintended BSs 



as noise. Hence, in this work, we assume SUD at the MS 
receivers. With SUD, it has been shown that transmit beam- 
forming is optimal for MISO IC in and (5). For the 
two-user case, Jorswieck et al. f6} proved that the Pareto- 
optimal beamforming vectors can be represented as linear 
combinations of the zero-forcing (ZF) and maximum-ratio 
transmission (MRT) beamformers. Previous studies j7l| and 
over MISO-IC beamforming usually assumed a central 
processing unit with global knowledge of all the downlink 
channels, which may not be feasible in practical systems. To 
make the result more implementable, our work focuses on 
multi-cell cooperative downlink beamforming, which involves 
distributed computations based on the local channel knowl- 
edge at each BS. Such decentralized multi-cell cooperative 
beamforming problems were previously studied in |£l based on 
the uplink-downlink duality to minimize the sum transmission 
power. Furthermore, a heuristic decentralized algorithm was 
developed in J4] for multi-cell cooperative downlink beam- 
forming based on the iterative updates of certain interference- 
temperature constraints across different pairs of BSs. 

It has been discussed in |4| that quantifying the Pareto 
optimal points in the achievable rate region over MISO IC may 
boil down to solving a sequence of convex feasibility problems 
after certain transformations, where the feasibility problems 
can be recast as second-order cone programming (SOCP) 
problems as shown in flOl . In this paper, we propose two 
algorithms to solve the resulting feasibility problem parallelly 
or distributively. In the first parallized beamforming algorithm 
based on alternating projections, we assume a computation- 
power limited centtalized processing unit such that part of 
the computation duties need to be parallely conducted in 
each individual BS. In the second beamforming algorithm, 
localized sequential optimizations across the BSs are per- 
formed iteratively, where the need for a central processing 
unit is eliminated. Convergence in norm for both algorithms is 
established. Besides, a set of feasibility decision rules is estab- 
lished to implement our algorithms for practical engineering 
applications. 

The rest of the paper is organized as follows. Section [TT] 
presents the MISO-IC model for multi-cell downlink beam- 
forming, defines the Pareto optimality, and reviews the rate 
profile approach, which transforms the whole problem into 
solving a sequence of SOCP feasibility problems. Section 



[HI] proposes the parallelized algorithm to parallelly solve the 
SOCP feasibility problem based on the method of alternating 
projections. Section [TV] presents the distributed beamforming 
algorithm to solve the SOCP feasibility problem based on 
the method of cyclic projections. Numerical examples are 
provided in Section [V] with conclusions in Section [VT] 

Notations: Bold face letters, e.g., x and X, denote vectors 
and matrices, respectively. J and denote the identity matrix 
and the all-zero matrix, respectively, with appropriate dimen- 
sions. diagiXx, . . . , X n ) defines a block diagonal matrix in 
which the diagonal elements are X±, . . . , X n . (•) and (-) 
respectively denote the transpose and the Hermitian of a matrix 
or a vector. W mxn and C mx ™ denote the space of m x n real 
matrices and the space ofmxn complex matrices respectively. 
1 1 a; 1 1 denotes the Euclidean norm of a complex vector x. All the 
log (■) functions are with base 2 by default. Re (•) and Im (•) 
denote the real part and imaginary part of a complex argument 
respectively, [ai; . . . ; a n ] defines a vector that stacks a\ . . . a n 
into one column. By default, all the vectors are column vectors. 

II. System Model and Preliminaries 
A. Signal Model 

We address downlink transmissions in a cellular network 
consisting of M cells, each having a multi-antenna BS to 
transmit an independent message to one active single-antenna 
MS. With the assumption that the same band is shared among 
all BSs for downlink transmissions, the system could be 
modeled as a Af-user MISO-IC. Specifically, we assume that 
each BS is equipped with K transmitting antennas, K > 1. 
With the assumption of single-user detection at each receiver, 
it has been shown in J4] and O that beamforming is optimal 
to maximize the rate region. Hence, the discrete-time baseband 
received signal of the active MS in the ith cell is given by 
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where uji £ C A denotes the beamforming vector at the ith BS; 
ha £ C K denotes the channel vector from the ith BS to its 
intended MS, while hji £ C K denotes the cross-link channel 
from the jth BS to the MS in the ith cell, % ^ j; Sj denotes the 
symbol transmitted by the ith BS; and Zi denotes the additive 
circular symmetric complex Gaussian (CSCG) noise at the 
ith receiver. It is assumed that Zi ~ CM (0, of) and Zi's are 
independent. 

We assume that the ith receiver only knows channel ha, 
and decodes its own messages by treating interferences from 
all other BSs as noise. With SUD, the achievable rate for the 
ith MS is thus given as 

Ri = log 1 + . " U I . ( 2 ) 



E 



h 



H, 



where the maximum transmission power is limited as 



where P, is the power constraint at the ith BS. 



B. Pareto Optimality 

We define the achievable rate region for the MISO-IC 
to be the collection of rate-tuples for all MSs that can be 
simultaneously achievable under a certain set of transmit- 
power constraints: 
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< r. t < Ri(u3x,. . .,um), 

1 = l,...,M 
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The upper-right boundary of this region is called the Pareto 
boundary, since it consists of rate-tuples at which it is im- 
possible to increase some user's rate without simultaneously 
decreasing the rate of at least one other users. To be more 
precise, the Pareto optimality of rate-tuple is defined as follows 
0. 

Definition 1: A rate-tuple (7*1, . . . , tm) is Pareto optimal if 
there is no other rate-tuple {fx, ■ ■ ■ , fjvf) with (fx,.. ., fu) 
> {rx,...,r M ) and (f 1; ...,f M ) ^ (r x , . . . , r M ), with the 
inequality being component-wise. 

In this paper, we are interested in searching the beamforming 
vectors for all BSs that lead to Parato optimal rate-tuples. 



C. Rate Profile Approach 

The rate profile approach ifTTI is an effective way to charac- 
terize the Pareto boundary of MISO-IC J4), where the key is 
that any rate tuple on the Pareto boundary can be obtained by 
solving the following optimization problem given a specified 
rate-profile vector, a = (ax, ■ ■ ■ , cum)' 
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where a. satisfies that a, > 0, 1 < i < M, and Ef=i a * = !■ 
Denote the optimal objective value of Problem © as R* sum . 
As shown in 0), R* sum ■ a corresponds to a particular Pareto 
optimal rate tuple. Hence, by exhausting all possible values 
for a, solving Problem (0 yields the whole Pareto boundary. 
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D. SOCP Feasibility Problem 

Directly solving Problem (O is usually difficult due to 
its non-convexity. However given the fact that the objective 
function is a single variable, we could adopt the bisection 
search algorithm to efficiently find R* sum as shown in H. 
Specifically, we could solve a sequence of the following 



feasibility problems each for a given ro'. 
max 
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i = l,2,...,M, 
Kf^j, j = l,2,...,M. (6) 

Therefore if the above problem is feasible for r , it follows 
that R* um > ro; otherwise, R* um < tq. Hence, a bisection 
search over R sum can be done. However, Problem © is still 
non-convex. 

As shown in iflOl . we can adjust the phase of u>j in (0 
to make h^uji real and non-negative without affecting the 
value of h|fu;j|. Hence, by denoting ft = e aiT ° — 1, i = 
1,2, ... , M, Problem (|6]l can be recast as 

max 

i = 1,2,..., M, 

0, i = l,2,...,M, 
ll^j-H < J = 1,2,...,M. (7) 

We further define x = [u)i; • • • ;u>m',0], n% = 
[0; 0; ... ; 0; a t ], Si = [ ■ ■ ■ I K • • • f with S,x = 
u> h i = 1,2, ...,M, and A l = diag (hg, fcg, . . . , h^, 0). 
For convenience, we add a term ft OJj) to both sides of 
the first constraint in Problem (|7]) as 

(1 + ft) (h*^) 2 > ft \Y^ =l \h%^\ 2 + o^j , (8) 

where i = 1, . . . , M. Finally, with our newly defined variables 
and coefficients, we recast Problem (O as 

max 

X 

S.t. v/ft \\AiX + m\\ < y/1 + Pi (/if SiX) , 

i = 1,2,...,M, 
p T a; = 0, 

ll^-xll <V^, 3 = !,••• ,Af, (9) 

where vector p is of the same dimension as x with all zero 
elements except for the last one being 1, such that the last 
element of x is guaranteed to be 0. 

Consequently, Problem |9) is a SOCP problem, which 
can be efficiently solved by numerical tools |fl~2). However, 
directly solving Problem (0 requires a centralized algorithm 
running at a control center, which may not be desired in 
certain engineering applications. Accordingly, there are usually 
two motivations for seeking distributed algorithms: one is 
to decompose the computations into multiple sub-programs 
such that the requirement for the central processing power is 
reduced; and the other is to localize computations such that no 



central control facility is required. In Section [TTU and Section 
IIVI we propose two algorithms based upon the above two 
motivations, respectively. 

III. Alternating Projections Based Distributed 
Beamforming 

In this section, in order to reduce the requirement on 
processing power at the control center, we develop a down- 
link beamforming algorithm, termed as alternating projections 
based distributed beamforming (APB), to solve Problem (O 
parallelly in M sub-problems. With our algorithm, the only 
processing power needed at the central unit is to calculate an 
average value over all the localized solutions from the M BSs. 
The algorithm is iterative, where parallel optimizations across 
BSs are performed at each round. The convergence issue of 
APB is also studied in this section. 

A. APB Algorithm 

At the initialization stage, the computation-limited cen- 
tralized unit is assigned with the values for M, K, and 
Pi, ... , Pm- Then the central unit broadcasts the information 
to all BSs with an arbitrary initial point xq G C +1 . It 
is assumed that the ith BS has the perfect knowledge of 
the channels from all BSs to the ith MS, i.e., all hij's. 
Furthermore, all BSs operate according to the same protocol 
described as follows. At the ?ith round, we denote the solution 
vector that the central unit broadcasts as x n -±. Then at the 
ith BS, the corresponding problem is expressed as 

min ||x-x„_i|| 

X 

S.t. y/fcWAiX + m\\ < y/l + Bi(hf[SiX), 

p T x = 0, 

\\S jX \\ < v/P~ j = 1,2,..., M, (10) 

where x n _i = jj Y^Li x n-i> with 4-1 denoting the 
optimal solution for Problem (JTOj of the (n — l)th round at 
the ith BS. A rough description of the algorithm is depicted 
in Fig. Q] 




Fig. 1. APB Scheme 



Remark 1: Note that if Problem ( TTOb is infeasible at the 
kth BS (k € {1, ...,M}), we can directly claim that the 



associated Problem (O is infeasible and quit APB. As such, 
from now on we only focus on the cases where Problem (fTOb 
is always feasible at each individual BS, and run APB to check 
when the overall problem in (O is feasible and when it is not. 

With a feasible Problem ( fTOb . we need the optimal solution 

a) 

Xn to satisfy all the transmitter power constraints and the ith 
receiver's SNR demand. The reason why we keep all M power 
constraints at each individual BS is for that fast convergence, 
which can be observed from simulations. Since all the Pj 
values are typically predetermined in cellular systems, no extra 
system overhead is needed. In the second-order cone constraint 
of Problem ([Tol l, directly using the term h$ SiX implies that 
Im (ft? Six) = and h% S t x > 0. 

B. Convergence Analysis 

Since APB is iterative, the convergence is an important issue 
to address. The convergence of APB is formally stated as 
follows. 

(i) 

Proposition 1: As n increases, the optimal solution x n 
for Problem (fTOb converges in norm to the limit x % when 
Problem © is either feasible or infeasible. Furthermore, the 
averaged solution x n also converges in norm to x* satisfying 
that jj 2~2i=i x l = x*. In particular, if Problem © is feasible, 
all x l 's coincide in the same point x that lies in the feasible 
set of Problem (O with x = x*. If Problem (O is infeasible, 
S*'s do not coincide in the same solution. 

Proof: 1) For the case of Problem (0 being feasible, we 
have the following proof. 

We first introduce the concept of finding the closest point 
to some given point in a closed convex set and alternating 
projections. 

In mathematics, a Hilbert space H is defined with the inner 
product (x, y) and the induced norm ||x|| = y (x, x ). If S is 
a nonempty closed convex set in H, Riesz iTPJI states that 
each x £ H has a unique best approximation (or nearest 
point) P s (x) in S. That is, \\x - P s (x)\\ < \\x - y\\ , Vy £ 
S\{Ps(x)}. The mapping P s : H -> S is called the 
projection onto S, i.e., finding the closest point to x in a 
closed nonempty convex set. In this paper, we use the general 
Euclidean inner product definitions (x,y) = x H y in the 
complex space and (x,y) = x T y in the real space. 

Definition 2: Suppose C\ and C2 are two closed nonempty 
convex sets in H with corresponding projections Pi and P%. 
Let C = C\ n C 2 and fix a starting point xq £ H. Then the 
sequence of alternating projections is generated by 

Xi = P 1 X ,X 2 = P2X 1} Xs = PiX 2 , 

xn = P2Xn-i,xn+i = P\Xn, ■ ■ ■ 

Let Fi denote the feasible set of Problem (fTOb at the ith 
BS, Fi ^ 0, and F = P|f=i Fi ^ ^ note that F is exactly 
the feasible set of Problem (|9). Thus solving Problem (TlOb at 
the ith BS can be viewed as finding the closest point to x n -i 
in a non-empty closed convex set Fi, i.e., the projection of 
x n -i onto Fi. Next we transform the variable defined over the 
complex Hilbert space to a double-dimensioned real Hilbert 
space such that we can use some existed results in alternating 



projections. We transform x g C N to x e R^ JV by letting 
x = [Re(x);lm(x)], where N' = KM + 1. Similarly, we 
map the complex set Fi to a double-dimensioned real set F[, 
and map x n -i to a double-dimensioned real vector x„_i = 
[Rc(x n _i); Im(sc n _i)]. We rewrite Problem ( fTOb as 
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P,-, j = l,...,M, 



(11) 



where 

At = 



Re (Aj) -Im^) 
Im(A0 Re (A,) 



,fc£=[Re(fc£) -Jm{h 



Si =[Si S t } , df = [ Im (h») Re (hg) ] [ Si S> 

n t = [m;0],p= \p;p] . 

From the constraints of Problem (fTTI) . we observe that the 
feasible set F( is the intersection of a collection of second- 
order cones, some subspaces, and some norm balls, which is 
nonempty closed and bounded. Next we show how to trans- 
form our algorithm into a problem of alternating projections. 
Let's define two product sets: 



T:F[xF^ x-.-Fl 



and 



U : {(a, a,..., a) : a E M 2Ar '| 



Meanwhile we define two new variables Xk,yk £ R 2Ar M as 



-(1) -(2) 



, Vk 



\Xk\X k \ 



\x k }. (12) 



Obviously, x k ET, y k £ U . By the results of Pierra in lfT4l . 
we have the following two lemmas: 

Lemma 1: Solving Problems ( fTOb for i = 1,...,M in 
parallel at the fcth round is equivalent to projecting vector 
yk-i onto the closed convex set T and obtaining x k . 

M 

Lemma 2: Computing jj ^ x ( £' is equivalent to project- 

i=l 

ing x k onto U and getting y k . 

Therefore, APB can be interpreted as alternating projections 
between T and U. Note that the idea of Alternating Projections 
was first proposed by von Neumann in |[T5ll , where only 
subspaces are assumed as the projection sets. Then many 
researchers extended this technique to more general scenarios 
|[T6l . ITF71 . For alternating projections between two non-empty 
closed convex sets C\ and C2, Cheney lfl6l proved that 
convergence in norm is always assured when either (a) one 
set is compact, or (b) one set is of finite dimension. Since 
set T is bounded and our underlying Hilbert space is of finite 
dimension, both conditions (a) and (b) are satisfied. Therefore, 
APB always leads to strong convergence, i.e., convergence in 



norm, due to the facts that the numbers of cells and antennas 
are always finite. As shown in ifTTl . all x i7 s will coincide into 
the same point x that lies in F. 

2) For the case of Problem (0 being infeasible, we have 
the following proof. 

With F = 0, the convergence of APB is still equivalent 
to the convergence of alternating projections between T and 
U, where Cheney's results in Ifl6l are applicable in this case. 
Thus, the convergence in norm is still valid for infeasible cases. 
Besides, it is easy to verify that jj $Z i=1 x l = x*. However, 
x l, s do not coincide into the same point. 

We now complete the proof for Proposition Q] □ 

C. Practical Feasibility Decision Rules 

With convergence in norm for APB established, we now 
need to establish some practical feasibility check rules to 
correctly terminate APB when it converges. 

From Lemma 1 and Lemma 2, we know that the feasibility 
of Problem (0 is totally determined by whether T and U 
intersect or not. By Proposition 1, if Problem (0 is feasible, 
all convergent solutions * coincide at a common 

point x* which belongs to F. In this case, all optimal values 
of Problems ( fTCTb converge to 0. On the other hand, if any 
of the optimal values of Problems ( fT0T > do not converge to 0, 
Problem (0 is infeasible. Based on the above discussions, we 
develop the following APB terminating procedures: 
Step 1: We set two threshold parameters e and £. The selection 
of e and £ affects the effectiveness of the algorithm. 
Step 2: Initialization: Let Vi, 1 < i < M, be the optimal 
value of Problem ( fTOb at the ith cell in the current computation 
round, v*, 1 < i < M, be the optimal value Problem (TTOb at 
the ith cell in the previous computation round, and flag[z], 1 < 
i < M be the flags for the M BSs. At the beginning, we set 
vi, . . . , vm and flag[l], . . . , flag[M] all zeros. 
Step 3: Repeat: For i = 1, . . . , M, the ith BS solves Problem 
( [Tol l and compares Vi against v*. If \v.i — v*\ > e, we refresh 
Vi : Vi = v* and proceed to Step 4; if \vi — v*\ < e, we 
compare v* with £: If v* > £, we claim that Problem (0 is 
infeasible and stop; otherwise, we mark this cell as flag[i] = 1 
and proceed to Step 4. 

Step 4: If flag[i] = 1 for all i = 1, . . . , M, we claim that the 
Problem (0 is feasible, then stop. Otherwise, return to Step 
3. 

Remark 2: Note that here we applied several approxima- 
tions in making the decisions. First, we claim that Problem 
([Tol l at the ith BS converges when \vi — v*\ < e. Thus, v* is 
considered as the limit of ith BS's optimal solution. Second, 
we set £ as the threshold dividing zero and non-zero values: 
If v* > £, we consider the limit non-zero, and vice versa. In 
simulations, we usually set both e and £ small with £ ^> e. For 
example, e = 0.002 and £ = 0.1 are chosen for the simulation 
results in Section M 

IV. Cyclic Projections Based Distributed 
Beamforming 

In this section, to localize computations such that no central 
control unit is required, we propose a decentralized algo- 



rithm that practically implements the multi-cell cooperative 
downlink beamforming. It is still assumed that the ith BS 
in the cellular network has the perfect knowledge of the 
channels from all BSs to the ith MS. Similar to APB, we 
decompose Problem (0 to M sub-problems and compute 
them at M BSs individually. In particular, the M problems 
are solved sequentially at each round, and the algorithm 
proceeds iteratively, which is termed as Cyclic Projections 
Based Distributed Beamforming (CPB). 

A. CPB Algorithm 

A certain cyclic update order among the M BSs needs to be 
determined at the initialization stage, where the 1st BS sends 
its solution to the 2nd, . . ., the (M — l)th BS sends its solution 
to the Mth BS, and the Mth BS sends its solution to the 1st, 
in a cyclic fashion. At the beginning, the M BSs should obtain 
the values for M, K, and Pi, ... , Pm- The algorithm starts 
from the 1st BS, after choosing an arbitrary initial point Xq, 
it solves the following problem 



mm 

xec KM + 1 



Xq 



s.t. 



ft||Aisc + ni|| < y/l + Piih&Six), 
0, 

< y/Pj, i = l,2,...,M, (13) 



T 
p X 

\\SjX 



where the optimal solution for the above problem is labelled as 
x± and sent to the 2nd BS. Then the other BSs begin to solve 
their own problems sequentially according to the predefined 
order. In particular, at the nth round the ith BS (i > 2) solves 
the following problem 



mm 

a;eC ifM + 1 



S.t. y/fcWAiX + m\\ < y/l + Pi(h%SiX), 



p T x = 0, 



< 



\Sjx\ 



Pj, j 1.2 M. 



(14) 



where x 



(i-i) 



is the solution sent over by the preceding BS, and 
Xn is used to denote the newly solved optimal solution. For 
simplicity, we refer the problem in (TPfl i as a cyclic subproblem. 
Such a scheme is illustrated in Fig. |2] 

Remark 3: Obviously, the constraints in Problem (TPfl i and 
Problem (TlOb are the same. Therefore, we have the similar 
discussions as in Remark [Tj We assume that all the cyclic 
subproblems are feasible when CPB is executed; otherwise, 
we directly claim that Problem (0 is infeasible. 

B. Convergence Analysis 

We first introduce the concept of cyclic projections. 

Definition 3: Suppose C\ , C2 , ■ ■ ■ , C r are closed convex 
sets in the Hilbert space H with C = Pi^Ci, and let Pj be 
the projection for Cj, i = 1,2, ... ,r. The operation of cyclic 
projections is an iterative process that can be described as 
follows. Start with any point x E H, and define the sequence 




BSm-j 



Fig. 2. CPB Scheme 

(x n ) (n= 1,2,...) by 
x = x,xi = Pi(x ), . . . , and x n = P nmo dr (x n -i) , (15) 

where Pk(-) is the projection operator to Cfe. 

In the literature, Bregman |fl8l showed that the above 
sequence generated by cyclic projections always converges 
weakly to some point Wc (x) £ C provided that C ^ 0, 
and Gubin |[T9l et al. provided a systematic study over general 
cyclic projections including the case of C — 0. Based on these 
results, we have the following proposition. 

Proposition 2: As n increases, the optimal solution x"n of 
the ith BS's cyclic subproblem converges in norm to a limit 
x % that lies in F{. Moreover, if Problem © is feasible, all x l, s 
coincide in a common point x* that lies in F. If Problem (O 
is infeasible, a; l 's do not coincide in the same solution. 

Proof: It is obvious that the optimal solutions for cyclic 
subproblems in (TPfl i form a sequence of cyclic projections. 
Since weak convergence is always guaranteed lfl8l . by the 
equivalence of weak convergence and convergence in norm in 
a finite dimensional space ll20"l . we obtain Proposition 2. □ 

Remark 4: Note that the convergence proof of CPB is 
more general than that of APB since alternating projections is 
actually a special case of cyclic projections where the number 
of projection sets is two. 

C. Practical Feasibility Decision Rules 

For CPB, the algorithm termination rules are similar to that 
of APB, which is skipped here. 

V. Simulation Results 

The performance of APB is first simulated. In the simula- 
tions, we set M = 3 and K = 4. We set the power constraints 
as 15, 18, and 21, respectively, for the three BSs. In Fig. [3] 
we demonstrates the convergence behavior as described in 
Proposition 1. The three curves correspond to the required 
SNR /3j's as 5, 10, and 20, respectively. We observe that their 
asymptotic behaviors are similar. In Fig. |4] with a feasible 



choice of = 10, i = 1, 2, 3, we draw how the achieved SNR 
values approach the target values over iterations. If Problem 
(0 is infeasible, for example, when setting target SNR as 
[50 40 60], the SNR evolution curves are given in Fig. [5] 
where we see that none of the target SNRs are satisfied. 

For the performance of CPB, the simulation setup is exactly 
the same as that for APB. In Fig. [6] with a feasible choice 
of Pi = 10, i = 1,2,3, we see how the achieved SNR 
values approach the target values with less iterations needed 
compared with Fig. |4] 

We observe from multiple simulations that the convergence 
speed of CPB is much faster than APB. 

VI. Concluding Remarks 

In this paper, based on alternating projections and cyclic 
projections, we have developed two optimal distributed beam- 
forming schemes to cooperatively solve the SOCP feasibility 
problem that is the key for quantifying the Pareto optimal 
points in the achievable rate region of MISO interference 
channels. The convergence in norm for both algorithms was 
established, which was further verified by numerical simula- 
tions. 
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Fig. 3. APB: || x n — x n — 1|| decreases 
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Fig. 4. APB: Achievable SNR tuple increases, on setting [10 10 10] 
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Fig. 5. APB: Achievable SNR tuple increases, on setting [50 40 60] 
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Fig. 6. CPB: Achievable SNR tuple increases, on setting [10 10 10] 
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